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We examine an effectively attractive quasi-one-dimensional Bose-Einstein condensate of atoms con- 
fined in a rotating toroidal trap, as the magnitude of the coupling constant and the rotational 
frequency are varied. Using both a variational mean-field approach, as well as a diagonalization 
technique, we identify the phase diagram between a uniform and a localized state and we describe 
the system in the two phases. 
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Effectively one-dimensional clouds of trapped atoms at 
very low temperatures are interesting systems because of 
a number of reasons. First of all, highly anisotropic traps 
allow us to realize this situation experimentally [1]. For 
example, in elongated traps, when the oscillator energy 
transversely to the long axis of the trap is much larger 
than the strength of the interaction, the degrees of free- 
dom along this direction are frozen out and the system is 
effectively one-dimensional. In addition, such systems are 
interesting from a theoretical point of view since they al- 
low us to check one-dimensional models which have been 
developed over the years. 

Recently Strecker et al. [2] and Khaykovich et al. [3] 
have managed to create and observe bright solitons - lo- 
calized blobs of atoms - in quasi-one-dimensional clouds 
of atoms after switching the effective coupling constant 
from positive (corresponding to an effective repulsive in- 
teraction between the atoms) to negative (corresponding 
to an effective attractive interaction) with use of the so- 
called Feshbach resonances [4]. References [5] have ex- 
amined these experiments theoretically. 

Inspired by these experiments, Refs. [6] and [7] have 
studied the behavior of a cloud of atoms confined in a 
toroidal trap. Their basic conslusion is that while the 
density of such a system is homogeneous for an effective 
repulsion between the atoms, when the interaction be- 
comes attractive, there is a critical value of the coupling 
constant below which the atoms form a localized blob. 
An easy way of understanding this effect is that the ki- 
netic energy that results from the Heisenberg uncertainty 
principle favors a homogeneous density, whereas an at- 
tractive interaction between the atoms favors the forma- 
tion of a localized density. 

Extending their study of Ref. [6], in a more recent 
paper Kanamoto et al. have examined the same problem 
when the toroidal trap is rotated [8]. While in a non- 
rotating torus there is only one parameter to vary, i.e, 
the ratio between the interaction and the kinetic energy 
(which we call 7) , in the present problem there is an extra 
degree of freedom, namely the rotational frequency of the 
torus f2. Using the mean-field approximation, as well as 
numerical diagonalization of the Hamiltonian, Kanamoto 
et al. have identified the phase diagram that separates 
the phase of uniform density from the phase of localized 



density in the f2 — 7 plane. 

In our study we attack the same problem using dif- 
ferent methods, and our results are in agreement with 
those of Ref. [8] . In the first method we use the mean- 
field approximation and calculate the order parameter 
variationally. In the second method we truncate appro- 
priately the Hamiltonian of the system and diagonalize 
it analytically using a Bogoliubov tranformation. Our 
results are also relevant with those of Lundh et al. [9], 
where the rotational properties of an attractive Bose gas 
that is confined in an anharmonic potential have been in- 
vestigated. As shown in this study, the effective potential 
felt by the atoms is toroidal and the angular momentum 
is carried cither by exciting the center of mass, or by cre- 
ating vortex states. One may argue that the localized 
phase of our problem corresponds to the center of mass 
excitation, while the uniform phase corresponds to the 
vortex phase [10] . 

Let us therefore consider a Bose-Einstein condensate of 
N atoms confined in a toroidal trap which has a radius 
R and a cross section S — irr 2 , where we assume that 
rCfi. Measuring the length in units of R and the energy 
in units of K 2 /2MR 2 , with M being the atom mass, the 
Hamiltonian H in the rotating frame of reference is 



H= alO 
Jo 

1 



$\0){L-nfj>(9) 



(i) 



Here ijj(0) is the field operator, L = —id/ 89 is the op- 
erator of the angular momentum, with 9 being the az- 
imuthal angle, 2f2 is the angular frequency of rotation, 
and Uo = 8TraR/S, with a being the scattering length for 
elastic atom-atom collisions. In the above expression the 
energy has been shifted by the term fl 2 for convenience. 

Let us start with the mean-field approximation. 
Within this scheme the system is described by the order 
parameter ^{9). Expanding ip(9) in the basis of plane- 
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wave states <pi(6) = e /v27r, we write 
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where we treat the coefficients c; as variational parame- 
ters. As shown in Refs. [11], the physical quantities of a 
system described by the Hamiltonian of Eq. (1) are given 
by the phase winding number J, 

J = [SI + 1/2], (3) 

where [x] is the largest integer that does not exceed x, 
and by the angular frequency relative to J, i.e., SI — J. 
Therefore the dominant term in the expansion of Eq. (2) 
is the one with I = J, where J is given by Eq. (3) for 
a given value of SI. When the density of the gas is uni- 
form, the occupancy of all other states is zero. However 
when the density is localized more states contribute to 
ip. Including the I — J ± 1 states we write, 

ip(9) = cj-x 4>j-\ + cj 4>j + cj+i <j> J+1 , (4) 

where in order for tp to be normalized, |c,/_ x | 2 + \cj\ 2 + 
|cj + i| 2 = 1. Close to the phase boundary the dominant 
terms in the order parameter are the ones with I — J,J± 
1, since other states (with I = J±2, etc.) have relatively 
higher kinetic energy and for this reason we can neglect 
them. As a variational method, our approach is more 
reliable close to the phase boundary, while further away 
from it, one needs to include more basis states. 

We proceed by expressing the energy per particle ej in 
the state of Eq. (4) in terms of the coefficients c; , which 
are then determined by minimizing e.j with respect to 
them [7,12]. We get that 

ej= (J-Sl) 2 \c.,\ 2 

+(J - 1 - ftflcj-il 2 + (J + 1 - 0) 2 |c. /+1 | 2 

+( 7 /2)[| cj _ 1 | 4 + m 4 + |4 +1 | 

+4| CJ _ 1 | 2 | CJ | 2 + 4| CJ | 2 | CJ+1 | 2 + 4| CJ _ 1 | 2 | CJ+1 | 2 

+2c 2 C }_ lC } +1 + 2(c 2 )* Cj _ iCj+1 ], (5) 

where 7 = NUq/(2-k). If we write ci — |c;|e* e ', there 
are three phases involved in the problem. One phase is 
arbitrary, the other one is free (reflecting the rotational 
invariance of the problem) and we choose it here so that 
the maximum of the density is at 9 = 0, and the third one 
is chosen so that the energy is minimum. To minimize 
ej we choose 8j-i + — 29 j = in order to make the 
last two terms in Eq. (5) as large as possible [7]. Using 
the normalization condition we eliminate cj thus getting 

ej - 7/2 -(SI- J) 2 = | CJ _!| 2 [1 - 2(J - SI)] 
+ \c J+1 \ 2 [l + 2(J-Sl)] 
+7[-|c j _ 1 | 4 -| Cj+1 | 4 + | Cj _ 1 | 2 + | Cj+1 | 2 
-|cj-i| 2 |c, 7+ i| 2 + 2|cj_i||cj+i| 
-2|c J _ 1 | 3 |c 7+1 | -2|c J _ 1 ||c J+1 | 3 ]. (6) 

When the density of the system is uniform, |c,/±i| = 0, 
and the energy per particle is ej ~ 7/2 + (SI — J) 2 , in 
agreement with the exact solution of the mean-field ap- 
proximation [8]. 



From Eq. (6) we now determine the phase boundary 
and the order parameter. Minimizing Eq. (6) we get two 
coupled equations, 

1 T 2(J - SI) + 7 [(A ±1 + 1)(1 - |c J±1 | 2 ) - 2\c.j Tl \ 2 

-3|c J+1 ||c J _ 1 |]=0, (7) 

where A = |cj + i|/|cj_i|. Close to the phase boundary 
\c,j±i \ is either zero, or < 1. This observation allows us 
to write Eq. (7) as 

l-2(J-fi)+ 7 (l + A) =0 
1 + 2(J - Q) + 7 (1 + 1/A) = 0, (8) 

and in order for a solution to exist, 7 = — l/2 + 2(J— SI) 2 , 
which is the equation that gives the phase boundary. Fig- 
ure 1 shows this curve and is periodic in SI, with a period 
of SI = 1. For 7 larger than this critical value the state 
is uniform, while for smaller values it is localized. Ref- 
erence [8] has derived the same result by examining the 
stability of the exact mean-field solution. Furthermore 
the transition between the two phases is of second order, 
since the order parameter is continuous across the phase 
boundary. Expanding we find that for 7 < —1/2 and 

\c l \ 2 = l -±^[l + a l (Sl-J)], (9) 
77 

where = —18/7 and aj +1 = 38/7, while |c/| 2 = 

i-|o-i| 2 -|o+i| 2 . 

Going even further, we have solved the two coupled 
Eqs. (7). Figure 2 shows the density of the gas \?p\ 2 for 
7 = —0.55 and —0.45 in each graph, and for \Sl — J\ = 
0, 1/4, and 1/2. For the smaller value of 7 the density 
is more peaked, as one expects intuitively. On the other 
hand the difference between the two states gets smaller 
as one gets further away from the phase boundary. Com- 
paring our variational solution with the exact one (within 
the mean-field approximation) for 7 = —0.55 of Rcf. [8], 
we see that in general the difference is small and it in- 
creases as \Sl — J\ increases. This is expected, however, 
since for higher values of \S1 — J\ one is further away from 
the phase boundary and more states have to be included 
in the expansion of the order parameter. Furthermore, 
when \Sl — J\ = the solution coincides with the one of 
a non- rotating cloud [6,7], apart from the phase factor 
e lJd , while for \Sl — J\ = 1/2 the order parameter has 
a node, since at this point the phase jumps by ir and J 
increases by one unit [8] . 

Let us now turn to the second approach we have de- 
veloped, namely the diagonalization of the Hamiltonian 
H. Following the same steps as before, we expand the 
field operator and for the same reasons as before, 

4>(9) = cj-i 4>j-i + cj 4>j + cj+i 4> J+ \, (10) 

where in this case q is the annihilation operator of an 
atom with angular momentum I. When the density of the 
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gas is uniform, the states I = J ± 1 are equally occupied 
and their occupancy goes to zero in the N — > oo limit. 
In addition the occupancy of other states (/ = J ± 2, 
etc.) goes more rapidly to zero (as N — > oo), since these 
states have relatively higher kinetic energy and for this 
reason we neglect them. Writing the Hamiltonian (1) in 
second-quantized form, 

H = ^(l - Sl) 2 cjci + ^ ^2 Ac\c m C n 5m+n-k-U (H) 
I klmn 

we notice that it can be diagonalized with a Bogoliubov 
transformation in the subspace of states with I — J, J±l, 
as we show below [12]. We work with the basis vectors 

\m) = \Nj_ u Nj,N J+l ) = \(J-ir,J N - 2m ,(J+l) m ), 

(12) 

where Ni is the occupancy of the state with angular mo- 
mentum I. The states |m) have been constructed to have 
N atoms, N — N, N J units of angular momentum, 
J2i l Ni = NJ, and an equal occupancy of the I = J ± 1 
states (i.e., m atoms in each of them). The diagonal 
matrix elements are, 

(m\H\m) = m[(J - 1 - ft) 2 + ( J + 1 - ft) 2 ] 

+ (N - 2m)(J - ft) 2 + - 2m)(N - 2m - 1) 

+2m(m - 1) + 8m(N - 2m) + Am 2 }. (13) 

When the system is on the side of the "uniform state" 
(in the limit N — > oo) then m is of order unity, which 
allows us to neglect the terms of order m 2 compared to 
those of order iV 2 and Nm (on the contrary, m becomes 
of order TV when a localized state forms), 

(m\H\m) «iV(J-ft) 2 +7(A r -l)/2 + 2m(l + 7). (14) 

We also get for the off-diagonal matrix elements, 

(m\H\m + 1) = ^2yJ(N - 2m)(N -2m- l)(m + l) 2 

47T 

»7(m + l) (15) 

in the limit m < N. From Eqs. (14) and (15) H can be 
written as 

H - j(N - l)/2 - iV(J - ft) 2 = 

(1+7X4-1^.7-1 +c\ +1 c J+1 ) 

+7(4-iCj+i + cj-icj+i). (16) 

The above Hamiltonian can be diagonalized with use of 
a Bogoliubov transformation [13], i.e., introducing the 
operators 

b = \iCj_j_ + A 2 cj+i and d = A 2 cj_i + Aic^ +1 . (17) 

Following the usual tricks, we get that the eigenvalues of 
H are given by 



£n b ,n d -j(N- l)/2- N(J-n) 2 = 

-( 7 +l) + v /TT2^(l+n 6 + n d ), (18) 

where rib, are the eigenvalues of the number operators 
b'b and cftd respectively. To get the above solution, 7 has 
to be larger than —1/2. For N — > 00, one can identify 
the term Eo^/N = 7/2 + ( J — ft) 2 as the total energy of 
the condensate with uniform density, in agreement with 
our variational approach and with the exact solution of 
the mean- field theory, Rcf. [8]. 

Furthermore, one can see that for the states |m), 
(b'b — cftd)\m) = 0, which implies that rib — rid — 
0,1,2,... Therefore, the first excited state of the sys- 
tem is (measured from the lowest energy) 2^/1 + 27 as 
Eq. (18) implies, if rib — n d = 1- This result is in agree- 
ment with Ref. [8]. 

Up to now the only condition we have found in order 
for a solution to exist is that 7 > —1/2. Turning to the 
stability of our solution against small perturbations, we 
consider two states where one atom with angular mo- 
mentum I = J — 1 has been transferred to the state with 
I = J + 1, or vice versa, 

|m ± ) = |(J - l) mTl , J N ~ 2m , (J + l) m±1 ). (19) 

The states \m^} have N atoms, and they carry an angular 
momentum that differs by two units compared to that of 
the state \m), J2i ^1 — NJ±2. One can follow the same 
steps as before and get that the energy is 

5JU d - 7(^V - l)/2 - N(J ft) 2 + 4(J - ft) = 

-(7 + l) + y/l + 2 7 (1 + n b + n d ). (20) 

In this case (b'b — d^d^m^) = ±2|m ± ) and thus the 
lowest energy of the system is that with rib — 2, rid = 0, 
or vice versa. In order for our solution to be stable, the 
lowest energy of £„ b n<j has to be higher than that of 

y/T+¥f>2\n-j\, (21) 

which can also be written as 7 > —1/2 + 2(J — ft) 2 , in 
agreement with the expression we derived earlier within 
our variational scheme. If Eq. (21) is not satisfied, even 
an infinitcsimally small anisotropy in the trapping po- 
tential suffices to convert the density of the system from 
uniform to localized. 

Finally we calculate the depletion of the condensate 
AN, 

AN = + = - 1, (22) 

which is singular only when J = ft. More precisely, close 
to the phase boundary and for small values of | J — ft|, 
AN w 1/(4| J - ft|), but it is finite otherwise. These 
results are also in agreement with those of the mean-field 
approach [8]. 
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To summarize, we examined a one-dimensional Bosc- 
Einstein condensate of atoms confined in a rotating torus, 
using two independent methods, namely a variational 
mean-field approach, as well as an exact diagonalization 
of the Hamiltonian. Based on these models, we examined 
the two phases which appear of uniform and localized 
density, and we identified the phase boundary which sep- 
arates them, as the strength of the interaction and the 
frequency of rotation vary. Our results are in agreement 
with those of the exact mean- field approach and with nu- 
merical diagonalization of the Hamiltonian of the system 
[8]. 
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FIG. 1. The phase boundary between uniform and local- 
ized density, 7 = -1/2 + 2(J - D.) 2 . The curve is periodic in 

with a period equal to 1. The stars denote the six cases 
considered in Fig. 2. 
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FIG. 2. The density of the trial order parameter |V>| 2 , 
Eq. (4) , calculated within the mean-field approximation, for 
|fi - J I = (top) \n-J\ = 1/4 (middle), and \Q - J\ = 1/2 
(bottom) with 7 = —0.55 in the solid curves and —0.45 in the 
dashed curves. 
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